Pair-correlation analysis of HD 10180 reveals a possible 
planetary orbit at about 0.92 AU 
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Since the observation in the first half of the nineties of exoplanets by precise timing measurements of a 
microsecond pulsar pQ as well as the observation of a Jupiter-like exoplanet around a main-sequence star 
f-^ ^ [2] around 500 exoplanets have been observed by techniques such as microlensing, astrometry, transiting, 

and imaging : 3j. Relatively recently, exoplanets with near Earth-like orbit have been reported [4] as well 
as planets with near Earth-like size [5]. The surge in exoplanet research and their observation has also 

OO given renewed impetus to discussions of planetary system formation [5] and to the century-old discussion 

CN of possible distance laws or regularities [HI 151 HO]. 

i— i Analyses for distance regularities akin to the Titius-Bode rule have suffered from what at times can be 

seen as arbitrary numbering schemes for planets (for a discussion of the undesirable in having to introduce 
additional parameters or other degrees of freedom, see Ref. [5]). We have suggested a pair-correlation 
analysis that does not involve any numbering schemes for the planets [TT]. Further, an analysis of the 
I Solar System demonstrated the correlation between the orbits of the planets extending far beyond nearest 

neighbour orbits. In this letter, we analyze the pair-correlations between the logarithmic positions of the 
5_i exoplanets in HD 10180. This exoplanetary system has six planets (plus one unconfirmed) surrounding a 

sunlike star [5]. The results of the analysis show that the six planets have a pair-correlation with nearly 
equally spaced peaks revealing the regularity of this system and its long-ranged correlation. Interestingly, 
the pair-correlation of the six planets reveals that a seventh orbit consistent with this group of six planets 
I exists. There is no indication for an eighth orbit. A modeling of the pair-correlations shows this additional 

>■ consistent orbit to be at 0.92±0.05 AU. 

The correlations between the positions of the exoplanets are determined from the pair-correlation 
function, which reads 

in 

ON 

o 
o 



P(A) = J p(x + A)p(x)dx. (0.1) 



Here p{x) is a distribution that measures the positions, Xi, of the planets on a logarithmic scale. For p 
we use a sum of Gaussian distributions so that p{x) = ^2 oii exp - ^ - ^) / 2<T , where Xi is the logarithmic 
position of planet i, Xi — log(aj/10 6 km), is the semi-major axis, and a = 0.075 corresponding to a 
• i-h full-width-at-half-maximum uj — 0.18. All weights, a*, are chosen as unity, a, = 1, as the pair-correlation 

does not depend much on the choice of o.i |11) . 

H 

The function P(A) is symmetric around the central self-correlation peak at A = 0. For the general 
case involving N planets there will be N(N — l)/2 peaks to either side of the self-correlation peak. With 
the Gaussian smearing above this will lead to broad wings. However, if the peaks are equidistantly spaced 
then the N(N — l)/2 peaks will coalesce into a smaller number of peaks. For the case of six equidistant 
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planets there will be five side-peaks indicating probable distances for the nearest neighbour, next nearest 
neighbour etc. 

In Fig. la is shown p(x) for HD 10180, i.e. the sum of the six corresponding Gaussian functions, and 
in Fig. Ib is shown the corresponding pair-correlation function, P(A). Fig. lb display a well structured 
set of nearly equidistant peaks. This is a manifestation of positional correlations in an exoplanctary 
system. Interestingly, there are six peaks to either side and not the expected five considering that the 
curve is generated using six planet positions. This suggests a seventh position. 

A first guess for the seventh position can be found from an inspection of the logarithmic positions 
of the exoplanets in Fig. la. Two candidate positions reveal themselves - between planets HD 10180 / 
and g, and between planets g and h, see Table 1. Fig. 2 corresponds to the first choice and Fig. 3 to 
the second. Fig. 2b shows enhanced correlations when compared to Fig. lb, while Fig. 3b shows a 
significant destruction of correlations when compared with Fig. lb. We therefore conclude that only a 
position between HD 10180 / and g is consistent with the six observed positions (c-h). By observing the 
calculated correlation pattern for a mesh of positions between HD 10180 / and g, the optimum position 
of the additional orbit is found to be at 0.92 ± 0.05 AU. 

Recently, an unconfirmed object close to the HD 10180 star has been reported [8]. When this object 
is included in the analysis (Fig. 4) the positional correlations are to a large extend destroyed, see Fig. 
4b. Therefore, this unconfirmed object is found not to belong to the same group of planets in the sense 
that its position is not in correlation with the six other reported planets. The seventh object is also 
significantly closer to the star (0.022 AU). 

The strength of the above analysis is that P(A) does not depend on any specific numbering of the 
planets. For example, it is independent of an interchange of the number index given to the planets. The 
calculated pair-correlation reveals a seventh possible orbit following a simple inspection of its number of 
peaks. When the origin of the type of orbit-orbit correlations are better understood we envision that 
it will be possible to extract quantitative numbers describing the planetary system under analysis and 
move beyond the qualitative inspections of pair-correlations function reported in this letter. 
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Table 1: Planets orbiting the sunlike star HD 10180. Here, ai is the semimajor axis, mass is in multiples of 
Jupiter mass, and Xi is ln(a/10 6 km). The data for ai and mass are from [8]. The two bottom rows are discussed 
in the text. 



Name 




a { 1 [10 6 km] 


Mass /[Mj] 


X i 


HD 10180 










HD 10180 b 


unconfirmed planet 


3.329 


>0.00425 


1.203 


HD 10180 c 


exoplanct 


9.589 


>0.04122 


2.261 


HD 10180 d 


exoplanet 


19.29 


>0.03697 


2.959 


HD 10180 e 


exoplanet 


40.38 


>0.07897 


3.698 


HD 10180 f 


exoplanet 


73.74 


>0.07520 


4.301 


HD 10180 g 


exoplanet 


212.7 


>0.06733 


5.360 


HD 10180 h 


exoplanet 


508.6 


>0.20262 


6.232 


HD 10180 f2 


consistent orbit 


138 




4.927 


HD 10180 g2 


inconsistent orbit 


329 




5.796 
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Figure 1: (A) Logarithmic position of six planets, not including the unconfirmed HD 10180 b. (B) 
Correlation function for these six planets. Notice that there are six side peaks rather than five. A linear 
regression of the positions of these six peaks gives an R 2 value of 0.996012. 




Figure 2: (A) Logarithmic position of the six orbits and a seventh at 138 x 10 6 km (f2) between / and g. 
(B) Correlation function for these seven orbits. A linear regression of the positions of the six peaks gives 
an R 2 value of 0.99861. 




*i A 
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Figure 3: (A) Logarithmic position of the six orbits and a seventh at 329 x 10 6 km (g2) between g and 
h. (B) Correlation function for these seven orbits. 
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